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Abstract 


The elastostatic plane problem of an infinite strip contain- 
ing two symmetrically located internal cracks perpendicular to 
the boundary is formulated in terms of a singular integral equa- 
tion with the derivative of the crack surface displacement as 
the density function. The solution of the problem is obtained 
for various crack geometries and for uniaxial tension applied to 
the strip away from the crack region. The limiting case of the 
edge cracks is then considered in some detail. The fundamental 
function of the integral equation is obtained and a numerical 
technique for solving the singular integral equations with this 
particular type of fundamental function which is characteristic 
of the edge cracks is described. The stress intensity factor 
for the complete range of net 1 i gament-to-wi dth ratio 0£a/h<J 
is calculated. The results also include the solution of the 
edge crack problem in an elastic half plane. 


1 . INTRODUCTION 

Within the past decade the plane elastostatic problem of an 
"edge crack" has been considered by many investigators using a 
variety of techniques. Among the notable contributions, we may 
mention [1,2] where the Wiener-Hopf technique is used to solve 
the problem, [3,4] where a variety of crack problems are formu- 
lated as integral equations in terms of some weight functions, 

[5] in which the problem is treated by using the mapping 
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technique, and [6,7] where the finite element methods are used. 

In this paper the problem is reduced to a singular integral 
equation with the derivative of the crack surface displacement 
as the density function. The method offers certain flexibilities 
which become very useful in dealing with the cracks in nonhomo- 
geneous materials. One of the primary aims of the paper is to 
show that the numerical techniques developed for solving singular 
integral equations with weight functions w(x) = ( 1 -x ) a ( 1 +x ) e , 

(-1 < Re(a,8) < 0) [8] may be extended to treat the problems in 
which formally a or 8 is zero, which invariably is the case for 
the edge cracks. 

First the problem of an infinite strip with two col linear 
cracks perpendicular to the sides will be considered. The limit- 
ing cases of this problem are a single crack when the distance 
between the inner crack tips becomes zero [9,10], and the edge 
crack when the distance from the outer crack tips to the sides 
becomes zero. The results for the edge crack in a half plane [2] 
will also be given largely for the purpose of verification of 
the technique. 

2. THE INTEGRAL EQUATIONS 

Consider the strip problem shown in Figure 1. Let the 
geometry of the medium and the external loads be symmetric with 
respect to x=0 and y=0 planes. Thus, using symmetry considera- 
tion and the superposition technique, the singular part of the 
solution may be obtained by solving the problem under the 
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following boundary conditions : 


u(0,y) = 0 , 

tf xy ( 0,y) = 0 , 

(0<y<~) , 

(1 .a,b) 

o xx ( h ,y) = 0 , 

o xy (h,y) = 0 , 

( 0 <y <°° ) , 

( 2 . a , b) 

a yy (x ’"> = 0 • 

Q 

X 

<< 

X 

8 

ii 

o 

v» 

( 0 <x£h ) , 

(3 . a ,b ) 

a xy (x,0) = 0 , 

( 0<x£ h ) , 


( 4 ) 

cfyyU, 0 ) = -a(x) 

, ( a<x<b) , 



§7 v(x ,0) = 0 , 

( 0<x<a , b<x£h ) , 


( 5 . a , b ) 

b a 

/ Iy v ( x , 0 ) dx = 0 

5 


(6) 


where u and v are the x and y-components of the displacement 
vector. The conditions (5.b) and (6) are clearly equivalent to 
v(x,0) = 0, (0<x<a, b<x<h). The solution of the problem satis- 
fying the necessary field equations and the condi tions ( 1 ) , (3) 
and (4) may be expressed as [8] 

u(x,y) = - | / (^- - py )e" py sin(px)dp 

0 H 

_ °o _ 

- f / {j[f(s) - !S y~ g(s}]sinh(sx) 

+ xg(s)cosh{sx) }cos (sy)ds , 


If the external load is the uniform tension a 0 as shown in 
Figures 1 and 2, then in (5a) a(x) = a 0 . However, in formulating 
both problems it is assumed that the crack surface traction may 
be a function of x. This may arise from symmetric but irregular 
loading of the strip and the half plane. 
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cos(px)dp 


V(x,y) - | 7 SM (H±l + py).-P» 

0 p 

+ I 7 {lc f (s) + *ir- g(s)]cosh(sx) 
o 5 * 

+ xg(s)sinh(sx)}sin(sy)ds , (7.a,b) 

where k = 3-4v for plane strain and k= (3-v)/(1+v) for general- 
ized plane stress, v being the Poisson's ratio. In (7) m, f, and 
g are unknown functions and may, in principle, be determined from 
the remaining boundary conditions (2) and (5). From (7) the 
stresses may be obtained as 

2 ^ a xx = " \ f m(p)(l-py)e' py cos(px)dp 

- — / [f ( s ) cosh(sx) + sxg(s)sinh(sx)]cos(sy)ds , 
o 

1 2 °° 

° yy = - £ f m(p) (l+py)e' py cos(px)dp 

+ I / {[f(s) + 2g ( s ) ] cosh ( sx ) 

o 

+ sxg(s)sinh(sx)}cos(sy)ds , 

1 O 00 

2y CT xy = " f / ypm(p)e' py sin(px)dp 

? 00 

+ £ f Uf(s) + g(s)]sinh(sx) 
o 

+ sxg(s)cosh(sx)}sin(sy)ds , (8.a-c) 

where y is the shear modulus. Substituting from (8) into (2) it 
is found that 
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f(s)cosh(sh) +shg(s)sinh(sh) = - ? cosh (ph) dp , 

* o (p 2 +s 2 ) 2 

[f(s) + g(s)]sinh(sh) + shg(s)cosh(sh) 

= ~ J - P 2 m ^^2 sin (P h ) d P • ( 9 . a , b) 

o {p +s ) 

The only unknown function is now m(p) which is determined 
from the mixed boundary conditions (5). (7.b), (8.b) and (9), 

substituted into (5), give a set of dual integral equations for 
m(p). However, a more direct approach to solve the problem would 
be defining a new unknown function G(x) as 

G( x )=-^v( x ,°), ( 0 < x < I h ) (10) 

and from (5.b) and ( 7 - b ) observing that 
+ i b 

- m(p) = / G(t)sin(pt)dt . (11) 

3 

Thus, substituting from (8.b), (9) and (11) into (5. a), we obtain 
an integral equation to determine G(x) of the following form: 
b 1 + 

/ H(x,t)G(t)dt = - Trcr 0 , (a<x<b) . (12) 

3 

It may easily be shown that at x=t the kernel H(x,t) has a Cauchy- 
type singularity. Separating this singularity, it is found that 

/ Ct^x + tTx + k ( Xst ) ' k(x,-t)]G(t)dt = - tto 0 , 

3 

(a<x<b ) , (13) 

where 

k(x,t) = / K(x,t,s)e"^ h_t ^ s ds , 
o 
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K(x,t,s) = e S ^{-[1 + (3+2sh)e 2s ^]cosh(sx) 

- 2sx sinh(sx)e _2sh - [2sx sinh(sx) 

+ (3 - 2sh + e" 2sh )cosh(sx)][l - 2s(h-t)]}/(l 

+ 4she _2sh - e' 4sh ) . (14. a, b) 

The index of the singular integral equation is +1; hence, its 
solution is determinate within an arbitrary constant which is 
determined from (6), or by (10), from 
b 

/ G(x)dx = 0 . (15) 

a 

Similarly, for a semi -i nf ini te plane shown in Figure 2 we 
find [11] 

/ G f-l dt + / N(x,t)G(t)dt = - tto 0 , (c<x<b) , (16) 

c c 

N(x,t) = ~ tTx + " 6x 2 “ ~~ 2 i > / G ( t )dt = 0 . (17. a, b) 

t+X (t+x) 2 (t+x) 3 c 

For 0<a<b<h and 0<c<b<°° the kernels k(x,t), k(x,-t), and 
N(x,t) are bounded, and hence (13) and (17) are ordinary singular 
integral equations with fundamental functions 

w(x) = [(x-a) (b-x)]" 1/2 , w(x) = [(x-c)(b-x)]" 1/2 (18. a, b) 

and may be solved in a straightforward manner (see, e.g., [8]). 

3. LIMITING CASE: THE EDGE CRACKS 

In the limiting case of a=0, b<h, the problem reduces to 
that of a single crack (-b,b) in an infinite strip which was 
considered in [9]. For this case, noting that G(x) = -G(-x), 
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(13) becomes 


/ I 777 + k(x,t)]G(t)dt = - ^7 , (-b<x<b) , (19) 

where the bounded function k(x,t) is given by (14). The funda- 
mental function of (19) is w(x)= (b 2 -x 2 ) ^ 2 and the solution 
subject to 
b 

/ G(x)dx = 0 (20) 

-b 

may be obtained by following any one of the standard techniques. 

In the other limiting case, i.e., for a>0, b=h, the problem 
becomes one of two edge cracks in a strip. In this case, through 
an asymptotic analysis, it can be shown that as both x and t 
approach the end point h the kernel k(x,t) in (13) becomes un- 
bounded and hence influences the singular nature of the solution. 
Since the integrand K(x,t,s) in (14) is bounded and continuous 
everywhere in 0<s<°°, the unbounded terms in k(x,t) will be the 
consequence of the asymptotic behavior of K(x,t,s) for s-*-«,. 

Thus, adding and subtracting the asymptotic value of K(x,t,s) to 
and from the integrand in (14. a) and using the relation [12] 

J s m e -s(2h-t) fsinh(sx)| ds = d^_ J e -s ( 2h- t ) fsi nh ( sx ) 1 ds 
o (cosh ( sx ) J dt m o (cosh ( sx ) J 

d m r 1 f x |i 

dt m L( 2 h - 1 ) 2 - x 2 ( 2 h - 1 )_ ’ (21 ^ 

we obtain 
h , 

/ [777 + k s (x>t) + k f (x,t)]G(t)dt = - a Q TT , 

(a<x<h) , (22) 
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k s (x,t) = [-1 + 6(h-x)j^ - 2(h-x) 2 ^] t . ( 2 h - x ) . 

00 . 

k f (x,t) = (t+x)' 1 - k(x,-t) + J[K(x,t,s) - K oo (x,t,s)]e"^ h_t ^ s ds , 

o 

K^CxjtjS) = e~ stl {-cosh ( sx) - [1 - 2s(h-t)][2sx sinh(sx) 

+ (3-2sh)cosh(sx)]} . (23.a-c) 


The first two terms of the kernel of (22) give a generalized 
Cauchy kernel. To obtain the fundamental function w(t) of the 
integral equation, let us first proceed in a routine manner [13, 
Chapter 4] and investigate the possibility of power singularities 
by defining 


( 24 . a-c ) 


G ( t ) = r(t)w(t) , w ( t ) = (h-t)“ a (t-a)" 3 , 

F(z) . i / M dt . 1 ( r(t)e 1,ia dt 

where r(t) is Holder-continuous in a<t£h and (t-h) a (t-a)^ is 
any definite branch which varies continuously in a<t<h. In (24) 
it is assumed that 

r(a) f 0 , r ( h ) f 0 , 0<Re(a,6)<l . (25) 

Noting that the point (2h-x) is outside the cut (a,h) on the real 
axis, (24) gives the following asymptotic relations [13]: 


F ( 2 ) 


= r(a)e 


TTi B 


1 


r(h) 


1 


(h-a) a sin7r6 (z-a)^ (h-a)^simra ( z - h ) a 0 


+ F„(z) , 


1 f G(t)dt _ r(a) cotTrg r(h) cotna „ r , x 

? l t G - ( (ih-x> * F < 2h - X > * - 77 r ie ) . TT^Ta + F 2 (*> 
a v ' (h-a) p smira (h-x) 
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(26. a-c) 



where the functions F.(z), (j=0,l,2), 

J 

a possible exception of the end points 
which [13] 

I F j ( z } I < ~~T P k ' P k ' Re(o.,e) 

Z-C k l 


are bounded everywhere with 
c k , (c-j =a , c 2 = h) near 


> ( j = 0,l ,2; k = l , 2 ) , (27) 


C k and p k being real constants. Substituting now from (26) into 
(22), multiplying both sides first by (x-a)^ and letting x=a, 
and then by (h-x) a and letting x=h, it is found that 


r(a)cot7r3 = 0 , 

j\ C2(l-a) 2 - 1 - cosira] = 0 . (28. a, b) 

(25) and (28. a) give 6 = 0.5 which is the well-known result. For 
r ( h ) f 0 the only possible root of (28. b) is a=l which is unac- 
ceptable. It is, therefore, clear that the function G(x) does 
not have a power singularity at the end point x=h. Before pro- 
ceeding with the determination of the fundamental function w(t), 
it should be pointed out that (28. b) is identical to the charac- 
teristic equation resulting from the asymptotic behavior of the 
solution for a 90-degree wedge. For example, if the wedge 
(0<r<°°, 0<6<tt/2) is subjected to boundary tractions 


a ee + ia re = 0 * 


(9 =-| , 0<r<°°) , 


CT 00 + ia r0 = » ( 8 = 0> r ! <r< r 2 ) , 

it can be shown that 

a,, a J r 2 i c + i°° . X-l 

TO IF u e< r *°> ' ’ ^ 2 FT — DTXT 


(29. a, b) 


D ( X ) = [2(1-X) 2 - 1 - cos7rX]/si nuX , 


(30. a, b) 
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where D(\) is the characteristic function. 


Since the function G(x) must have an integrable singularity 
at x=a, and since the foregoing analysis indicates that it cannot 
have a power singularity at x=h, the next step is to look for a 
solution of the following form: 


G(x) = r 0 (x)(t-a)’ e , (0< Re ( 8 ) < 1) (31) 

where r Q (a) and r Q ( h ) are bounded and nonzero. Defining again 
the following sectionally holomorphic function 


F(z) 



r Q (t)dt 
(t-z) (t-a) B 


(32) 


we obtain 


r n (a)e 
F(z) = -2 


TT i 8 


•1"* (zV + l09(2 - h) + F3<z) ' 


(x-a) 


tt ( h - a ) 

rAh) 


i / r«n ' F ' 2h -*> ■ l09(h - , ° + Fs<x> 


7 A / - F 6(x) 


G(t)dt 


(t, -^ )2 \ • F 7 (»> • < 33 — > 
where the functions F 3 and F^ are bounded everywhere except 
possibly at the end point z=a where they have a behavior similar 
to (27); in particular they tend to definite limits as z h and 
x + h. The functions Fg, Fg and Fy are bounded and tend to defi- 
nite limits as x -► a and x+h. Now substituting from (33) into 
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(22) it is seen that the logarithmic singularities are cancelled 
and the power singularity (x-a) can be eliminated by requiring 
that cotuB = 0, which again gives 3=1/2. We may thus conclude 
that the solution of (22) is of the form (31). Formally then 
the weight function of the integral equation is (h-t)°(t-a)~^ 2 
with the related polynomials P n ^*"^ 2 ^(t). Similarly in (16) 
for c = 0 by letting G( t) = r-j ( t) (b-t)"^ (where r-j is a bounded 
function), it may easily be shown that 3 = 0.5. It should be 
pointed out that for c=0 and b=h the general solution of (16) 
and (22) no longer contains an arbitrary constant; therefore the 
conditions (15) and (17. b) are not necessary for a unique solu- 
tion which is also suggested by the physics of the problem. 

4. ON THE NUMERICAL SOLUTION OF THE INTEGRAL EQUATION 

Consider the singular integral equation (16) which, for 
c =0, has a generalized Cauchy kernel. To solve the integral 
equation numerically by techniques similar to that described in 
[8], we first normalize the interval (0,b) by defining 

x = t/b , £ = x/b , 

G(t) = <Ht) - r,(t)(b-t)- 1/2 = R(t)(1-t)- 1/2 . (34) 

Using (34) (16) may be expressed as 

/ M(C,t)0(x) dx = - ffcr(b£) , (0<C<1) . 

o M 

M(£,t) = (t- 5)" 1 + bN ( b£ , bx) . (35. a, b) 

To solve (35) numerically the most practical technique appears 
to be to extend the definition of <j>(x) in an appropriate manner 
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into the interval (-1,0) and to use the corresponding Jacobian 

integration formula with the related orthogonal polynomials 
(-1/2 -1/2) 

P n 5 ; (t). Noting that $(0) / 0, an appropriate extension 

of <j> may be an even continuation as follows: 

4>(t) = p(t)(1-t 2 )’ 1/2 , p ( t ) = p ( -t ) , ( — 1 < t< 1 ) * 

p(t)(1+t)' 1/2 = R ( x ) , ( 0 < t< 1 ) . (36) 

With (36) (35) becomes 

} MU,T)*(T)dT = \ } M(5,|t|) TTO(bO , 

0 2 -1 ( 1 - T 2 ) 1/2 

(0<?<1). (37) 

Using now the integration formula corresponding to the 

weight function (l-x 2 )" 1/2 [8] and from (13), (16) and (17) 

observing that x=0 is a zero of M(£,x), we obtain 

1+1/ . •, 2n + l 

- - 2,1, A iP( T i)M(C j .|T i l) 

n 

= I A 1 p(x i )M(C i ,T i ) , ( j = 1 n ) , (38) 

i =1 1 1 J 1 

T 2n+1 ^ T i ^ 0 » T i 00 s ^ * A i = 2n+T 

U 2n^p 0 * cos ^2n+T^ * ( 39 ) 

(38) provides a system of n linear algebraic equations to deter- 
mine the unknowns p(x.j)> (0<x.<l, i = l,...,n). Once the density 
function G(t) is obtained the stress intensity factors and the 
crack surface displacement may be obtained from (see Figure 1) 

k (a ) = lim /2(a-t) a (x,0) = lim /2(t-a) G(t) , 
t->a t->-a K 
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5. NUMERICAL RESULTS 

As a first example consider the half plane problem shown 
in Figure 2. The numerical results given in this section are 
obtained for a uniform pressure a(x)=a Q applied to the crack 
surface which corresponds to the uniform tension applied to the 
strip (Figure 1) or the half plane (Figure 2) perpendicular to 
and away from the plane of the crack. 


Table 1. Stress intensity factor at x = b for the 
half plane (a Q =(b-c)/2, c Q =(b+c)/2). 


c /a 
0 0 

1.05 

1.01 

1.001 

1.0 

1 .0 (Ref .2) 

■HQK 

1 . 2540 

1.3303 

1.3987 

1.5869 

1 . 5861 

o n /a~ 
0 0 


The results for the half plane are given in Table 1 which also 
shows the result given in [2]. The agreement appears to be quite 
good. Further results for larger values of c Q /a 0 and the values 
of k(c) may be found in [11]. Figure 3 shows the crack surface 
displacement for the edge crack, c Q /a 0 =1 and for an internal 
crack, c /a =1.01 obtained from 

d 

v(x,0) = - / G ( t ) dt . (41 ) 

x 
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The normalizing constant v Q shown in the figure is 

v 0 = o 0 a 0 (l+«)/4y . (42) 

The results for the strip are given in Tables 2 and 3, and 
in Figures 4 and 5. Table 2 gives the results for two symmetri- 
cally located collinear internal cracks shown in Figure 1. In 
limit when a=0 the results reduce to that found in [10] for a 
single (central) crack. As the size of the net ligaments 2a and 
(h-b) approach zero, the corresponding stress intensity factors 
go to infinity. These limiting values are indicated in Table 2 
by an arrow. From the Table it is seen that the rupture of one 
of these ligaments causes a sharp increase in the stress inten- 
sity factor at the other crack tip (for example, compare the 
lines 1,2 and 5,6). 

Table 2. Stress intensity factors for collinear internal 
cracks in a strip (Figure 1, a =(b-a)/2). 


a/h 

b/h 

k(a) 

Mb) 

0 

0.4 

(+ 00 ) 

1.5690 

0.1 

0.5 

1.1746 

1.1169 

0.2 

0.6 

1.1102 

1.0961 

0.4 

0.8 

1.0984 

1.1250 

0.5 

0.9 

1.1290 

1.2278 

0.6 

1.0 

1.6080 

(■* 00 ) 

0 

0.8 

(-> 00 ) 

2.5680 

0.1 

0.9 

1.6730 

1 .7451 

0.2 

1 .0 

2.1769 

(+ 00 ) 

0.5 

0.95 

1.1960 

1.4711 

0.5 

0.98 

1.2713 

1.9008 

0.5 

1.0 

1.6228 

(-*- 00 ) 
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Table 3. Stress intensity factor for the edge cracks in 
a strip of finite width (b = h, a 0 = (h-a)/2, 
net section stress: a-| = a Q h/a). 


a/h 

k(a) 

°o^o 

Ma) 

o-j /a" 

0 

-> 00 

|= 0.63662 

0.1 

2.9467 

0.62510 

0.2 

2.1769 

0.61572 

0.3 

1 .8744 

0.60738 

0.4 

1.7136 

0.59361 

0.5 

1.6328 

0.57728 

0.6 

1.6080 

0.55703 

0.7 

1 .5970 

0.51749 

0.8 

1 .5915 

0.45014 

0.9 

1.5883 

0.33690 

1.0 

1.5869 

0 


Table 3 gives the results for the edge cracks which are also 
shown in Figure 4. The limiting values of these results are 
that of the edge crack in a half plane for a/h -> 1 (see Table 1) 
and infinity for a/h->-0. For the case of (h/a)=°°the closed 
form solution is given by (e.g., [14]) 


%y( x >°) = - 7== 


( -a<x<a ) , 


k (a ) = lim /2(a-x) a (x,0) = 


x+a 


yy 


■n/a 


(43 . a , b ) 


where P is the resultant load (per unit thickness) acting along 
the y-axis. Defining now a net section stress o-j by 

(44) 


CT 1 ' 2a = a o a 


P _ _ h 
a 

for a/h + 0, the stress intensity factor becomes 
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0.63662 . 


(45) 


k(a) = 2 = 
a ] /a 77 

The two stress intensity factor ratios shown in Table 3 and 
Figure 4 are related by 


Kill „ iiii_ (1 . 


a-| /a 


cr /a7 
o o 


(46) 


Note that for h/a+°°, for a fixed a k(a) is proportional to 
a Q h or a Q a 0 . Hence the ratio k(a)/(a 0 /a^)-*- °° as h/a-*-°° or 


a Q /a 


Examples of the crack surface displacement v(x,0), (a<x<b) 
for the internal and edge cracks in the strip are shown in Figure 
5. The normalization constant v Q used in Figure 5 is also given 
by (42). 
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